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Abstract. The magnification bias of large-scale structures, combined with galaxy biasing, leads to a cross-correlation of distant 
quasars with foreground galaxies on angular scales of the order of arc minutes and larger. The amplitude and angular shape 
of the cross-correlation function wqg contain information on cosmological parameters and the galaxy bias factor. While the 
existence of this cross-correlation has firmly been established, existing data did not allow an accurate measurement of wqg yet, 
but wide area surveys like the Sloan Digital Sky Survey now provide an ideal database for measuring it. However, wqg depends 
on several cosmological parameters and the galaxy bias factor. We study in detail the sensitivity of wqg to these parameters 
and develop a strategy for using the data. We show that the parameter space can be reduced to the bias factor b, Qq and erg, and 
compute the accuracy with which these parameters can be deduced from SDSS data. Under reasonable assumptions, it should 
be possible to reach relative accuracies of the order of 5%-15% for b, fio, and erg. This method is complementary to other 
weak-lensing analyses based on cosmic shear. 
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1. Introduction 

Weak gravitational lensing by large-scale structures gives rise 
to a magnification bias on distant quasar samples. Quasars be- 
hind matter overdensities are magnified and therefore prefer- 
entially included into flux-limited samples. Area magnification 
dilutes quasars on the sky and counter-acts the magnification to 
some degree, but the number count function of bright quasars is 
sufficiently steep to provide a huge reservoir of faint quasars to 
be magnified above the flux limit. Bright quasars in flux-limited 
samples thus occur preferentially behind matter overdensities. 

Galaxies are biased tracers of dark matter. If the bias is 
positive, their number density is higher in matter overdensi- 
ties. Since these overdensities act as gravitational lenses for 
background sources, they give rise to purely lensing-induced 
cross-correlations between background quasars and foreground 
galaxies. 

Typical angular scales for this effect range from a few arc 
minutes to one degree. The lensing effect of individual galaxies 
does not matter on such scales. Group- or cluster-sized haloes 
contribute to the correlation only on scales below about one 
arc minute. On larger scales, the expected signal is caused 
by the large-scale matter distribution only. Measurements of 
the lensing-induced QSO-galaxy correlation function therefore 
have the potential to constrain the dark-matter power spectrum, 
several cosmological parameters, and the galaxy bias. 



Many such measurements have been undertaken 
(e.g. Bemtez & Martinez-Gonzalez 1995, Norman & Impey 
1999, Bemtez et al. 2001, Norman & Impey 2001). The exis- 
tence of highly significant large-scale QSO-galaxy correlations 
has been convincingly demonstrated (cf. Fugmann 1990; 
Bartelmann & Schneider 1993, 1994; Bartelmann et al. 1994; 
Rodrigues- Williams & Hogan 1994; Benitez & Martinez- 
Gonzalez 1995; Norman & Impey 1999), but the measurement 
of the amplitude and shape of the cross-correlation function is 
still highly uncertain (see Bartelmann & Schneider 2001 for 
a review). This is mainly because homogeneous, sufficiently 
deep galaxy surveys of large portions of the sky have so far 
been unavailable. 

Upcoming wide-field surveys, above all the Sloan Digital 
Sky Survey (SDSS; York et al. 2000), will provide huge, homo- 
geneous samples of quasars and galaxies covering a substantial 
fraction of the sky. The detection and analyses of QSO-galaxy 
cross-correlations with high signal-to-noise ratio is thus com- 
ing within reach. We investigate in this paper which parameters 
can best be constrained with such a measurement, and what ac- 
curacy we can expect to achieve. 

The main route to extract cosmological information from 
weak lensing has so far been the exploitation of cosmic shear. 
This approach measures the gravitational tidal field of the dark 
matter between background galaxies and the observer. Despite 
the many difficulties in measuring exact shape parameters of 
faint galaxies, impressive results have been obtained in the 
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past year, demonstrating the power of the method (Bacon et 
al. 2000, Hoekstra et al. 2001, Kaiser et al. 2000, Van Waerbeke 
et al. 2001, Wittman et al. 2000). However, it depends upon the 
crucial assumption that galaxy ellipticities are intrinsically un- 
corrected (see Heavens et al. 2000 and references therein). 

Quasar-galaxy cross-correlations induced by weak lensing 
provide a complementary method for analysing weak lensing 
by large-scale structures which does not depend on this as- 
sumption, and does not need any shape measurement. It there- 
fore allows an independent and welcome cross-check of the 
cosmic shear results and, in addition, it allows the determina- 
tion of the galaxy bias factor. 

The ratio flo/b can be constrained using the ratio 
wqq/wqq, where w GG is the angular auto-correlation func- 
tion of the foreground galaxies (Bemtez & Sanz 2000). 
However, this method neglects all information contained in the 
amplitude and angular shape of wq G . Motivated by recent cos- 
mic shear results, we investigate here how well cosmological 
parameters can be constrained if the full information provided 
by quasar-galaxy cross-correlations is used. 

In this paper, we briefly outline the lensing theory of QSO- 
galaxy cross-correlations and investigate on which parameters 
the cross-correlation function depends. We show that the di- 
mensionality of the parameter space can be substantially re- 
duced. We then use simulated measurements to investigate how 
accurately the remaining parameters can be constrained. 

2. Theory 

The theory of the QSO-galaxy cross-correlation function 
caused by weak lensing was introduced and developed in sev- 
eral earlier studies. We can therefore be brief here and limit the 
discussion to the issues required later on. 

Let tlq and uq be the number densities on the sky of 
quasars and galaxies, respectively. The QSO-galaxy cross cor- 
relation function is then 



{[gQ(g) - n Q }[n G (8 + 4>) - n G }) 



nQ n G 



(1) 



where the average extends over all positions <fi and all direc- 
tions of 6. The average number densities of quasars and galax- 
ies are uq and n G , respectively. If there is no overlap in redshift 
between the quasar and galaxy populations, this correlation is 
exclusively due to lensing, and can be written (Bartelmann, 
1995) : 



w QG (0) = 2 (a -1)6(0) 1^(0) 



(2) 



provided that the magnification \i is weak, \fi— 1| <C 1. Here, 
a is the logarithmic slope of the cumulative quasar number- 
count function near the flux limit and w K s(0) is the cross- 
correlation function between the lensing convergence k and 
a suitably weighted projected density contrast 5. The conver- 
gence itself is a weighted projection of the density contrast, so 
that the correlation function w K s (0) is straightforwardly related 
to the projected dark-matter power spectrum. 

Equation ^ provides an operational definition of a mean 
galaxy bias factor b(9) on a given angular scale 0. As defined, 
it averages over all galaxies considered in the cross correlation 



at an angular separation from the nearest quasar, and thus 
involves averaging over redshift and possibly also morpholog- 
ical type, unless galaxies are selected according to their mor- 
phology. Being linear in the galaxy number density, the QSO- 
galaxy cross-correlation is insensitive to a possible stochastic- 
ity of the bias. In practice, measurements will determine the 
wqg within an angular range R around 0. We will abbreviate 
the average bias factor within that interval as 6r. 

Of course, the relation of the bias factor bji to commonly 
used definitions of the galaxy bias factor depends on how the 
bias depends on physical scale, redshift, and galaxy type. We 
will further discuss this issue in Sect. 2.2 below. 

2. 1. The k-6 cross correlation 

We now describe briefly how the cross-correlation function 
w K s between lensing convergence and projected density con- 
trast can be evaluated. More detail can be found in Bartelmann 
(1995), Dolag & Bartelmann (1997), Sanz et al. (1997) and 
Bartelmann & Schneider (2001). 

Since both k and 6 are weighted line-of-sight projections 
of the density contrast, we aim at an expression relating w K s to 
the matter power spectrum. 

Specifically, we can write k as the projection 



K (6) 



dwp K (w)S[f K (w)6,w] , 



(3) 



where w is the radial comoving distance along the line-of-sight, 
fx(w) is the comoving angular diameter distance and p K (w) is 
the projector 

Pk(w) = 77 ( — 
2 \ c 

I — Wq(w) j. , „ . (4) 



a(w) 



f K {w>) 



The ratio between the angular diameter distances represents 
the usual effective lensing distance, Wq(w) is the normalised 
distance distribution of the sources, in our case the quasars, and 
a(w) is the cosmological scale factor. 

Similarly, the projected density contrast 8 can be written as 



6(6) 



dw pg(w) 5[f K (w)9,w] 



(5) 



where ps(w) is the normalised distance distribution of the 
galaxies that are cross-correlated with the quasars. 

The correlation function w K s can now be related to the 
statistics of 6, 

w kS (9) = (k(c/>) 8(4 + 6)) 

= J dwp K (w) J dw' ' ps(w') 

x (8[f K (w)6,w]5[f K (w')(<f> + d),w']) . 

Using Limber's equation for the statistics of projected homo- 
geneous Gaussian random fields, inserting the Fourier trans- 
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form of the density contrast, and introducing the power spec- 
trum Pg(k), we find 

(a \ f A Pk(w)ps(w) 

x r^ Ps {hky w ) Mseh (6) 

The quasar-galaxy cross-correlation function, wqq(9) = 
2 (a— 1) b(9) w K g(9), describes the statistical excess of quasars 
around galaxies with respect to a Poisson distribution, or the 
excess of galaxies around QSOs, at an angular distance 9. 

2.2. The bias factor 

The relation of the average bias factor bu to the usual galaxy 
bias factor depends on the properties of biasing and of the 
galaxies considered. If biasing is linear and possibly stochas- 
tic, the bias factor as defined here averages over rb val , where 
6 var is the ratio of variances in the galaxy and dark-matter dis- 
tributions, and r is the linear correlation coefficient defined by 
Dekel & Lahav (1999). For nonlinear biasing, b^ averages over 
the biasing parameter b, which is the slope of the linear re- 
gression of the galaxy fluctuations on the density contrast and 
whose relation to other conventional measures of galaxy bias- 
ing becomes more involved (Dekel & Lahav 1999). In other 
words, the physical meaning of b^ is defined by Eq. (^), but its 
relation to other measurements of galaxy biasing depends on 
the biasing model and remains to be specified if such a relation 
needs to be established. 

Generally, the galaxy bias can depend on time, scale, 
and galaxy parameters like luminosity or morphological type. 
Numerical simulations using semi-analytic galaxy-formation 
models predict how the bias of different morphological galaxy 
types evolves with redshift (cf. Kauffmann et al. 1999), and 
such prescriptions can be used for establishing the relation be- 
tween bn and other biasing parameters. The impact of a pos- 
sible scale dependence of the bias on quasar-galaxy cross- 
correlations can also be taken into account. This was studied 
by Guimaraes et al. (2001), who used ratios between measured 
and theoretical galaxy power spectra for estimating the mean 
bias parameter in Fourier space, b(k). Naturally, such methods 
strongly depend on the bias model. 

For our application, the redshift dependence of the mean 
bias b can probably be neglected because the galaxies ex- 
pected to be cross-correlated with background quasars will be 
in a relatively narrow redshift range peaking near z ps 0.3. 
Furthermore, the scale dependence can be constrained. Bemtez 
& Sanz (2000) showed that the ratio wqq/wqq effectively 
measures J7o b/bl ar , without the need to know the cosmolog- 
ical parameters. For linear and deterministic biasing, 6 var re- 
duces to b and we have : 

wgg{8) b{9) 

for quasars and galaxies located at single redshifts zq and zq- 
Van Waerbeke (1998) extended this calculation to realistic red- 
shift distributions, for which the constant Q becomes a known 



scale-dependent function Q(9), but remains only weakly sen- 
sitive to cosmological parameters. Therefore, the method pro- 
posed by Bemtez & Sanz allows constraints on the scale de- 
pendence of the bias. 

We investigate here the possible constraints on cosmologi- 
cal parameters that can be derived from the angular shape and 
amplitude of wqq. For doing so, we assume that the bias is 
independent of scale on the angular scales we are considering, 
i.e. between 1' and 1°. A scale-dependent bias factor could be 
taken into account assuming that the variation with scale can 
be expressed by a smooth function f(9), so that 

b(0)=b R f(9) (8) 

Models or measured constraints on f(9) can then be included 
into the calculation of the cross-correlation function wqg(9). 

2.3. Expected galaxy overdensity 

We now evaluate Eq. (||) under the following assumptions. 
First, we choose a CDM power spectrum and use the formal- 
ism by Peacock & Dodds (1996) to approximate its non-linear 
evolution. This is crucial because the linear approximation un- 
derestimates the correlation amplitude on arc-minute scales 
by about an order of magnitude (Dolag & Bartelmann 1997). 
Unless specified otherwise, we use the cluster normalisation 
constraint cr 8 = 0.52 f)-°- 52 +°- 13n ° from Eke et al ( 19% ) 

Since we shall later need to compute the correlation func- 
tion many times to explore the parameter space, we assume 
that the redshift distributions for source quasars and foreground 
galaxies are delta functions, which speeds up the computation, 
but affects the amplitude of w K s(9) only by at most 20% com- 
pared to realistic distance distributions. We place all galaxies at 
redshift z — 0.3, and all quasars at redshift z = 1.5. 

Finally, we need a model for the quasar number counts as 
a function of flux. In the quasar catalogue of the Sloan Early 
Data Release (Schneider et al. 2001), the cumulative number- 
count function for quasars brighter than 19 th magnitude is well 
approximated by a power law with a logarithmic slope of a ~ 
2. 

Some illustrative results for wqq (9) are presented in Fig. |l} 
The curves in the left panel show the expected cross-correlation 
functions for nonlinearly evolving density perturbations in a 
flat, low-density universe (Slo = 0.3, ft\ = 0.7, cluster nor- 
malised with er 8 = 0.93; solid line), an Einstein-de-Sitter uni- 
verse (cluster normalised with er 8 = 0.52, dashed line; and 
normalised to er 8 = 1, dashed-dotted line). 

For comparison, the dotted curve was computed for an 
f?o = 0.3, J7a = 0.7, cluster normalised cosmology, but as- 
suming linear growth of the density perturbations. At small an- 
gular scales (below ~ 10'), most of the correlation amplitude 
is contributed by mildly nonlinear matter perturbations. The 
right panel shows the correlation function out to large angular 
scales, where it drops below zero. The break towards two de- 
grees occurs when the first zero of the Bessel function in Eq. (||) 
moves across the physical scale of the peak in the power spec- 
trum. Hence, the correlation signal is expected to drop sharply 
at these large angular scales. 
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Fig. 1. Left panel: The QSO-galaxy cross-correlation function wqq(6) is plotted for different cluster normalised cosmologies. A 
low-density flat universe with fio = 0.3, 57a = 0.7 is plotted in solid line; an Einstein-de Sitter universe is plotted in dashed line, 
and in dashed-dotted line for a universe normalised with erg = 1. Both curves were computed using the non-linearly evolving 
matter power spectrum. For comparison, the dotted line shows how the correlation function if we consider linear growth of the 
power spectrum for a 57o = 0.3, 57 a = 0.7 cosmology. Here we have assumed a = 2, a bias factor of b = 1 and a shape 
parameter of V — 0.21. The right panel shows how the correlation function drops at large scales before reaching negative values 
in a low-density flat universe. 



Figure |l] demonstrates that the quasar-galaxy cross- 
correlation function varies substantially if the cosmological 
model is changed. Therefore, it can be used as a tool to mea- 
sure cosmological parameters. We explore this issue in the next 
section. 

3. Dependence on cosmological parameters 

3. 1. Reduction of the parameter space 

Equation (^) illustrates how the quasar-galaxy cross-correlation 
function depends on various assumptions and parameters. The 
redshift distributions of the source quasars and the foreground 
galaxies enter through the functions p K (w) and ps(w), respec- 
tively. These distributions can be observed. Likewise, the log- 
arithmic slope a of the quasar number-count function can di- 
rectly be obtained through observations. 

Therefore, the remaining unknowns of the cross-correlation 
function are the dark-matter power spectrum and its growth, the 
cosmological parameters, and the bias factor of the galaxies 
relative to the dark matter. We now explore how sensitively the 
cross-correlation function depends on these inputs. 

First, we shall assume that the matter content in the 
Universe is dominated by cold dark matter. For a Universe con- 
taining pure cold dark matter, the location of the peak of the 
power spectrum is defined by V = 57o h (Bardeen et al. 1986), 
and it changes to T = 57o h exp[— 57b (1 + \2h / 57o)] in pres- 
ence of baryons (Sugiyama 1995). For reasonably low values 
of the baryon density parameter, the main effect of the preced- 
ing modification is merely a weak change of the first relation 



between T and 57 . For 57 b h 2 = 0.02 (Wang et al. 2001), we 
find that T ~ 57o h — 0.03. This detailed, weak dependence on 
the baryon density does not qualitatively affect the following 
results. We shall see in Sect. 4.2 how the parameter estimation 
changes if we ignore any relation between V and the other pa- 
rameters. 

The remaining free parameters appearing in the quasar- 
galaxy cross-correlation function are the cosmological param- 
eters 57o, 51a and h, the normalisation of the power spectrum, 
which is expressed by its variance as on scales of 8Mpc/ft,, 
and the bias factor bu. 



3.2. Independence on VL\ 

Interestingly, w K s(9) depends only very weakly on the cosmo- 
logical constant 57 a on scales larger than one arc minute. This is 
demonstrated by Fig. ||for different angular scales and different 
cosmologies. Above one arc minute, the contours tend to align 
parallel to the 57a axis. At an angular scale of 15 arc minutes 
for instance, the variation of w K s with 57a for a given density 
parameter 57 is less than 10%. For comparison, increasing 57 
from 0.1 to 1 changes the correlation amplitude by roughly a 
factor of four. Thus, focusing on intermediate and large angular 
scales, the correlation function can be considered insensitive to 
57a- This safely allows the reduction of the effective parameter 
space to 57 , erg, h and &R. 
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9 = 2 arcmin 0=15 arcrnin = 40 arcmin 




Fig. 2. The dependence of w K s on Hq and £7 a is shown for two different angular scales (6 = 2', 15' and 40')- Here we have used 
a cluster normalised cosmology, with r = Hq h and 6 = 1. On scales greater than one arc minute, the 57a -dependence is so weak 
that this parameter can safely be eliminated. 



3.3. Dependence on fl anda 8 

Regarding the dependence of w K g on fio and <jg, a rough in- 
spection of Eq. would give w K s °c £lo<J%, where the factor 
Slo comes from Poisson's equation, and the scaling oc erf re- 
flects the normalisation of the power spectrum. However, closer 
examination reveals a more complicated dependence. Due to 
the nonlinear evolution of the power spectrum, as does not only 
change the correlation amplitude, but also its shape. Similarly, 
fig affects the nonlinear evolution of the power spectrum, and 
therefore also changes the shape of w K g. 

Figure |^ illustrates the dependence of w K s on f^o and as, 
at three angular scales, for fixed values of h and 6r. For mea- 
surements taken at a single angular scale only, the dependence 
on £lo and as is degenerate. At small angular scales, for in- 
stance, an increase in can be compensated by a decrease 
in <7g (cf. the left panel of Fig. [|). However, this degeneracy 
can be broken if measurements at different angular scales are 
combined. At large angular scales, for example, an increase in 
fig would have to be compensated by an increase in erg (cf. the 
right panel of Fig. ^J). Therefore, the shape of the correlation 
function contains cosmological information, and it is thus im- 
portant to measure the correlation signal across a wide range 
of angular scales. The dashed line in Fig. ^| shows the cluster 
normalisation constraint given by Eke et al. (1996). 

We note that there is an interesting angular scale around 40' 
where the quasar-galaxy cross-correlation function becomes al- 
most independent of S7o- The dependence of w K s on Slrj and as 
can be understood as follows: 

- From Eq. (§, we see that the correlation function is roughly 
proportional to Orj x J sds P(s) Jq(s6). At small angular 
scales, the Bessel function is close to unity over a large 
range of wave numbers s, and the result of the integra- 
tion over s is roughly proportional to a\. Hence, we have 
w oc f2rj °8' whatever the value of T. 

- At large angular scales, the oscillations of the Bessel func- 
tion largely cancel the power at high wavenumbers, and 
thus the s integration is only sensitive to the low-s part 
of the power spectrum. Since this rising part of the power 
spectrum decreases rapidly with increasing V (or equiva- 



lently increasing Qq) an d keeping the normalisation con- 
straint fixed, the correlation amplitude decreases as 51o in- 
creases. 

- At intermediate angular scales, the various effects of chang- 
ing i7o cancel each other almost exactly, which leads to the 
insensitivity of the correlation amplitude to Hq at angular 
scales around 40'. 

The overall dependence of w K g on as and f2o reflects the 
physically motivated proportionality T oc f2rj> which is the sim- 
plest possible relation between these parameters. More general 
cases will be taken into account in Sect. 4.3. 

Our main conclusions from this section are that the depen- 
dence of w K s on Ha is so weak that it can be neglected, and 
the degeneracy between as and Hq can be broken by measur- 
ing w K s over a wide range of angular scales between a few arc 
minutes and one degree. Assuming T = Ho h, w K $ becomes 
almost independent of flo for angular scales near 40'. At that 
scale, the amplitude of the correlation function is proportional 
to 6<t|, and a measurement of this product is in principle pos- 
sible. We will address this point in Sect. 4.2. 

4. Expected measurements 

We now evaluate how accurately the remaining cosmologi- 
cal parameters can be constrained through measurements of 
quasar-galaxy cross-correlation functions. For that purpose, we 
simulate observations of quasar-galaxy cross-correlations ex- 
pected to be possible with the Sloan Digital Sky Survey data. 

4.1. Signal-to-noise estimate 

Since our goal is to constrain cosmological parameters from 
quasar-galaxy cross-correlations, we now assess a strategy for 
optimally using the data. Within narrow, ring-shaped bins of 
radius 9 centred on A^q background QSOs, the total deviation 
of the galaxy count from its mean is 

{N G!i -N G )^N Q N G wq_ g (6) , (9) 

i=l 
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Fig. 3. The dependence of w K s on Slo and as is shown for three different angular scales (9 = 2', 40' and 60' from left to right) 
and assuming r = f2o h. We kept h = 0.72 and 6 = 1 fixed for the plot. As the dependence on as and f2rj changes with angular 
scale, the degeneracy between these two parameters can be broken. At ~ 40', the correlation function is almost insensitive to Qq. 
The dashed line shows the cluster normalisation from Eke et al. (1996). 



where Nq is the mean galaxy count obtained from randomly 
selected fields, and Ng,{ the galaxy count in the ring centred 
on the i-th QSO. In a stack of such narrow, ring-shaped bins, 
the signal-to-noise ratio of the detection will be 



N Q N G w QG {0) 



(10) 



As we saw before, the quasar-galaxy cross-correlation function 
wqg (9) is larger than the cross correlation w K s between con- 
vergence and projected density contrast by a factor 2 (a — 1) b. 
The value of the bias factor 6r is expected from simulations 
fall near unity. However, the value of (a — 1) depends on the 
quasar sample. 

The measured cumulative quasar number-count function 
strongly depends on the selection criteria of the quasars. For 
our measurement, we first need a cut at low redshift to exclude 
any overlap between the QSOs and the galaxies. Then, for a 
given magnitude range, we need a clean and complete sample. 
Since most of the SDSS quasars are photometric, this implies 
cutting the sample in multicolour space for optimising the se- 
lection criteria. We leave this level of detail aside in this paper. 
We use the fitting function provided by Pei (1995) for estimat- 
ing the cumulative quasar number-counts as a function of mag- 
nitude, i.e. a broken power law with slopes of a ~ 0.64 for 
QSOs fainter than B ~ 19, and a ~ 2.52 for brighter QSOs. 

Since faint QSOs have a < 1, they are negatively-biased 
by the lensing magnification. Thus, taking all observed QSOs 
straightforwardly into account would tend to reduce or even re- 
verse the correlation signal. However, we can separate the de- 
pendences of wqq (9) on angular scale (which contains the cos- 
mological information) and on quasar magnitude and thus keep 
all information required to stack the correlation signal from all 
quasars in a sample. 

The Cauchy-Schwarz inequality provides optimal weights 
for quasars with respect to magnitude. We find that the signal- 
to-noise ratio is optimised by 



opt 



/ dm (a(m) - l) 2 N Q (m) 
f dm\a(m) — 1| yj Nq (m) 



where A*Q(m) is the number of quasars in the magnitude bin 
[m, m + dm], and a(m) is the corresponding slope of the num- 
ber counts. Therefore, we can consider the following effective 
correlation function 



w c s(0) = 2 & R (a off - 1) w K s(9) 



(11) 



where a e g is an effective slope describing the entire quasar 
sample, 



Qfefi = 1 



1 /dra(a(m) - l) 2 N Q (jn) 
/JVq J" dm \a(m) - 1| yjN Q {m) 



(12) 



The signal-to-noise ratio of the optimally weighted data can 
now simply be calculated using w e s instead of wqg in 
Eq. (|ic|). The effective correlation function w c s(9) is related 
to an average absolute deviation of the galaxy number counts 
in the vicinity of quasars from their mean, weighted appropri- 
ately. In effect, each magnitude bin gives a constructive signal 
to the cross-correlation, and the weighting ensures a maximised 
signal-to-noise ratio. To get an idea of a reasonable value for 
a e g, we use Pei's (1995) broken power-law fit for a sample 
of QSOs ranging in magnitude from m = 16 to ra = 21. 
Evaluating the integral in Eq. (p"l|), we find that a e g = 1.35. 

Some quasars have a magnitude in for which \a(m) — 1| « 
0, so they will only contribute noise to the final correlation 
function. Those quasars are down-weighted by Eq. (jll|), but 
the selection can be further restricted by keeping only quasars 
whose signal-to-noise ratio for the lensing effect exceeds a 
specified threshold. The sample is then constrained by the con- 
dition 



| a (to) - 1| > 



y/NQNG2b R W KS (9) 



(13) 



/ N G 2b R w KS (9) , 



The number of quasars satisfying this criterion is expected to 
be a large fraction of the total available sample. We thus neglect 
the sample restriction in the following, but it should be taken 
into account in any measurement of wqg from real data. 

In the previous section, we emphasised the importance of 
measuring the cross-correlation across a wide range of angular 
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scales with a good signal-to-noise ratio. Thus, the size of the 
ring-shaped bins around quasars in which galaxies are counted 
must change with the ring radius 0. At large radii, wide bins are 
needed, and the earlier narrow-bin calculation no longer holds. 
The general result for an annular bin is 

where Ue is a top-hat filter of radius 0. This integration over a 
disk replaces the zeroth-order Bessel function Jq(sO) in Eq. (^) 
by 27r6>Ji (s6)/s, and thus 

p K {w)p s (w) 



d 2 0W Kl5 (0) Ug((p) = I dw r .. 



Jk( w ) 



(15) 



x / dsP s 



w 0Ji(s6 



Jk{w) ' 

In order to adapt this model to a forthcoming application to 
the SDSS data sample, we will use the following parameters: 

- The galaxy density is t%q ~ 3 arcmin (Scranton et 
al. 2001) for SDSS galaxies down to r' = 22. It is shown 
that masking of stars or bad seeing is important at that mag- 
nitude limit, and the usable area of the survey is reduced by 
a factor 2/3. 

- The angular coverage is taken from 6 = 0.5' to 1° and 
separated in 15 logarithmically spaced bins. Below 0.5', the 
shot noise becomes important due to the fairly low galaxy 
number density. Above one degree, the amplitude of the 
correlation function drops rapidly as shown in Fig. |], and 
so does the signal of the galaxy overdensity. 

- The effective logarithmic slope of the quasar number-count 
function is set to a c s = 1.35. 

- The quasar number Aq is essential for determining the 
signal-to-noise ratio. Its value is continuously increasing 
with the size of the survey and will reach its maximum 
in 2005. The SDSS Early Data Release contains approx- 
imately 4000 spectroscopic quasars and 9000 quasar can- 
didates based on photometry alone. The expected quasar 
number of the final survey is on the order of 10 5 for spectro- 
scopically identified quasars, and can reach on the order of 
10 6 for photometrically selected quasars (G. Richards, pri- 
vate communication). Furthermore, we need a quasar sam- 
ple which is well separated in redshift from the galaxy sam- 
ple. For this purpose, we will only use quasars with z > 1. 
In order to account for these numbers and for the reduction 
of the usable survey area due to masking, we will consider 
an effective number of quasars Aq between 5 x 10 3 and 
5 x 10 5 . 

Photometric redshifts will be available for all quasars. The 
amplitude of wqq is mainly sensitive to the mean redshift 
of the considered quasar sample rather than its exact red- 
shift distribution. Therefore we can safely neglect any un- 
certainties in the redshift estimates. 



4.2. Constraints on 6(40') a| 

We saw earlier that assuming Y = fio h (or similar relations) 
causes the quasar-galaxy cross-correlation to be only weakly 



sensitive to Qq near the angular scale of ~ 40' (correspond- 
ing to an effective scale of roughly 8 h^, 1 Mpc), and this prop- 
erty can in principle be used to further reduce the parameter 
space. At angular scales near 40', the density perturbations 
probed by the correlation function are well within the linear 
regime, where the amplitude of the power spectrum simply 
scales oc cr|. At that scale, the correlation amplitude scales 
oc 6(40') erf to good accuracy, almost independently of the den- 
sity parameter fio and the cosmological constant SI a. Given 
Qrj> the shape parameter of the power spectrum is fixed by the 
Hubble constant. The product 6(40') a\ can therefore be deter- 
mined through the correlation amplitude only up to the level of 
uncertainty reflecting the uncertainty of the Hubble constant, 
and by the remaining weak dependence on fio- Figure Q shows 
how the amplitude of the quasar-galaxy cross-correlation func- 
tion near 40' changes with the Hubble constant and the density 
parameter fio. The overall uncertainty is on the order of 20%. 



8.e-4 



6.e-4 



4.e-4 



2.e-4 



O.e-4 



■0.4 



0.6 



0.8 



Fig. 4. w K $(£l ; 9 = 40') is plotted for different values of Y. 
The solid line is for T = fio h with h = 0.72, and the two 
dashed line are for h = 0.64 and h = 0.80, which are the 1-cr 
confidence limits given by the Hubble Key Project. Here we 
have used erg = 1. 

The variation of w K s with f2o can be reduced by narrow- 
ing the angular range from which the data are taken, however 
at the expense of the signal-to-noise ratio. Within an annulus 
with inner radius 30' and outer radius 50', where the ampli- 
tude of wqg(^o) changes by less than about 20%, the expected 
signal-to-noise ratio can be determined from Eq. ( |l4| ) with the 
appropriate top-hat filters, 



0r(5O' 2 - 30' 2 ) 



9 



10 5 



Nq\ 1/2 (ug\ 1/2 f 2(a c ff-l)fr(40') 



V 3 J 



0.7 



(16) 



where we have used the number of spectroscopic quasars ex- 
pected to be provided by the SDSS. Using photometric quasars, 
this signal-to-noise ratio might be improved. Therefore, the ac- 
curacy of a measurement of 6(40') u\ through quasar-galaxy 
cross-correlations is limited by the uncertainties in the Hubble 
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constant and the density parameter if f2rj is treated as unknown. 
We emphasise the fact that the only external information we 
use for this method is the Hubble constant. Any additional con- 
straint on 57 can substantially increase the precision of the 
6(40') erf measurement. 

4.3. Parameter estimates: Q , 08 andb R 

We now investigate possible parameter estimates using the in- 
formation coming from all angular scales. Using Eqs. ( 14) and 
( 15 ), we can simulate measurements of wqq(9) and then com- 
pute the uncertainty of the parameter constraints we wish to 
extract from the data. 

We saw above that the shape of the correlation function 
contains important information on the cosmological parame- 
ters. For achieving similar signal-to-noise ratios in all annular 
bins of the measurement, we use logarithmically spaced annuli 
around the quasars. In order to estimate the uncertainty of our 
parameters, it is worth noticing that 10 5 quasars spread evenly 
across the whole survey area are on average separated by about 
20'. Given our 15 logarithmically spaced bins covering the an- 
gular range introduced above, this means that the errors in the 
three outer bins would be effectively correlated. The others are 
only very weakly affected by possible correlations. Given the 
low amplitude of the cross-correlation, the Poisson noise domi- 
nates in our case. Hence, a simple x 2 of the difference between 
the simulated data and the model enables us to reasonably es- 
timate the uncertainty of the parameters. However, it has to be 
kept in mind that we are ignoring any angular variation in the 
bias, which might contribute to the parameter uncertainties. 

However, due to the large parameter space and some partial 
parameter degeneracies, the \ 2 minimisation yields allowed 
ranges for the complete parameter set which are fairly wide 
compared to the existing parameter constraints. An attempt to 
find constraints on the full parameter space (f2 , ct§, 51a, h, 6r) 
is therefore not promising. We have to include additional infor- 
mation for better constraining our parameter set. 

The prime advantage of the quasar-galaxy cross- 
correlations compared to other methods for estimating 
the cosmological parameters is its direct sensitivity to the 
galaxy bias factor. In addition, it is a weak-lensing based 
method which can directly measure £7o and provide an in- 
dependent cross-check of results obtained from weak-shear 
observations, as it does not rely on the assumption of intrin- 
sically uncorrelated galaxy ellipticities. The most accurate 
measurements of fio will be provided by CMB experiments. 
They rely on information coming from photons emitted at the 
recombination epoch (z ~ 1200) and depend on the curvature 
of photon trajectories as they propagate through the Universe. 
An independent measurement of Hq using only low-redshift 
information probes the consistency of the cosmological model. 

For the following, we will use two external pieces of infor- 
mation. The first is an estimate of the power spectrum normal- 
isation erg. For this, we use the cluster normalisation relation: 
ct 8 = 0.52r^ - 52+0 1300 (Eke et al. 1996), since there is no 
precise and direct measurement of erg. The second is an esti- 
mate of the peak location of the power spectrum T. We will use 



either a direct measurement, provided by galaxy surveys, or a 
measurement of the Hubble constant h, if we adopt the relation 
r = Slo h. Finally, the remaining unknown parameters of the 
problem are the density parameter f2o, the density parameter 
corresponding to the cosmological constant and the bias 
parameter 6r , of which £7 a can be ignored. 

Other reductions of the parameter space are of course pos- 
sible, but this one is particularly attractive because it allows a 
measurement of the galaxy bias and its dependences on scale, 
galaxy morphology, magnitude etc. In what follows, our first 
fiducial model has r = firj h, h = 0.72, in agreement with 
recent constraints from the Hubble Key Project (Freedman et 
al. 2001), and is cluster normalised. We will later drop this re- 
lation between V and Qrj- Hence, we compute 

2m m V- [SNu(O i -,n ,b R )-8N oba (e i -,Q ,bR.o)] 2 nn . 
X (fin, 6) = 2^ 2 

i a i 

where 5Nyi(9i) and SN ^ s (9i) are the overdensities of galax- 
ies in bin number i expected in the model and detected in the 
simulated observation, respectively, and <jj is the Poisson noise 
in that bin. 

The results are shown in Fig. |5| For this example, we have 
assumed a cosmological model with flo = 0.3, 51a = 0.7 and 
a galaxy bias factor of unity. The contours show the 68.3%, 
95.5%, and 99.7% confidence levels. Results for two different 
quasar sample sizes are presented, namely 10 4 and 10 5 quasars, 
from left to right. We see that a large quasar sample is required 
in order to at least confirm a low Hq. 

Previous quasar surveys listed only a few thousand ob- 
jects, therefore we expect the SDSS to provide us with the 
first quasar-galaxy cross-correlation measurement with a rea- 
sonable signal-to-noise ratio that allows the extraction of some 
cosmological information. The two contour plots show a near- 
degeneracy between £7 and 6r, for low fl . Indeed, by lower- 
ing Sl the correlation function increases due to the cluster nor- 
malisation. For low n , mainly the amplitude of the correlation 
function changes rather than its shape. Low values of the bias 
parameter can compensate for this change and still allow a good 
fit. Thus, assuming the relation T = Qrj h, quasar-galaxy corre- 
lations can provide an estimate of fio, but only poorly constrain 
6r. Figure || also shows that 6 can be measured with good accu- 
racy if the cosmological parameters are assumed to be known; 
A6 R « 15% for Q = 0.3. 

As we discussed before, it is possible to probe other pa- 
rameters with the quasar-galaxy cross-correlation function. For 
example, Fig. ^ shows what confidence intervals we can expect 
for 51o an d 08> assuming h — 0.72 and 6r = 1. For instance, 
an independent measurement of the bias can be obtained com- 
bining cosmic shear and galaxy counts provided similar galaxy 
samples are considered (cf. Schneider 1998 and van Waerbeke 
1998 for the theory, and Hoekstra et al. 2001b, in preparation, 
for a measurement using a joint analyses of the Red Cluster 
Survey and the VIRMOS-DESCART Survey). Using external 
information on 6r, we see that interesting constraints on 51o 
and as can be derived: we find A57o ~ 10% and Acts ~ 5%. 

The validity of the assumed relation T = Qq h (or more 
generalised relations taking baryons and neutrinos into ac- 
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Fig. 5. Expected constraints on the parameters CIq and b for cluster normalised universes. The input parameters used to create the 
simulated observation are (f^o, ^a, eg, T) = (0.3, 0.7, 1.0, 0.21), and the bias was set to unity. The panels show the 1-, 2- and 
3-cr confidence levels for quasar sample sizes of 10 4 , 10 5 , using the relation r = CIq h. The symbol "+" shows the location of 
the input model, and the symbol "x" shows the minimum of the x 2 - 




Fig. 6. Expected constraints on cosmological parameters. The contours show the 68%, 90%, and 95% confidence levels. The 
two panels show results for quasar sample sizes of 10 4 , 10 5 . The parameters of the input model are (Oo, ^a, cs, T) = 
(0.3,0.7, 1.0, 0.21), and b R = 1. 



count) depends on the matter and radiation content of the 
Universe, which is presently uncertain. For a more conserva- 
tive approach, one can drop this relation and use measure- 
ments for T. The two plots in Figure [7| show the correspond- 
ing changes of the confidence levels. In the ficr^R plane, we 
see that we then get interesting constraints of the parameters. 
As we can expect in the flo-as plane, the contours are not sig- 
nificantly affected, and an interesting sensitivity on Hq and eg 
remains. Similar contour plots can be obtained from cosmic- 
shear measurements (e.g. Van Waerbeke et al. 2001) by assum- 
ing a fixed value for T. As shown in Figure 0, quasar-galaxy 



cross-correlations measured with the SDSS will soon enable us 
to probe these parameters equally well as cosmic shear does 
with current deep imaging surveys. 

Regarding f^o and 6r, we saw from the different contour 
plots that we achieve much higher sensitivity by using external 
information on T. Table 1 summarises the uncertainties on Clo 
and 6r we can achieve, depending on the quasar sample size. 
We can reasonably expect the final usable quasar number to be 
at least of the order of 10 5 , which means an independent mea- 
surement of Slo wim an accuracy of ~ 17%, and of 6r with an 
accuracy of ~ 10%, should be feasible. 
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Fig. 7. Expected constraints on cosmological parameters. The contours show the 68%, 90%, and 95% confidence levels, for a 
quasar sample size of 10 5 , considering the fixed value T — 0.21. The parameters of the input model are (firj; ^A, 08, T) = 
(0.3, 0.7, 1.0, 0.21), and Br = 1. The frame in the left panel shows how the 1-er contour moves by considering V — 0.21, but 
having r trU e = 0.19 or 0.23. 



Clearly, the final goal of using quasar-galaxy correlations is 
probing the galaxy bias, but as these plots show, this method 
can also be used for constraining cosmological parameters. 



effective quasar number 


Afi A6 R 


1 x 10 4 
5 x 10 4 
1 x 10 5 


0.14 0.25 
0.07 0.12 
0.05 0.09 



Table 1. 1-er uncertainties expected for ilo an d bji as functions 
of the effective quasar number, using SDSS parameters and the 
input values ilo = 0.3 and Br = 1. Here, we consider a cluster 
normalised ACDM cosmology and a fixed value for V. 

The additional parameter V is now known to approximately 
20% (Dodelson et al. 2001), and its accuracy is expected to 
greatly improve with SDSS in the near future. However, this 
uncertainty adds to the uncertainty of flo and 6r. To see the 
implications, the frame in the left panel of Fig. ^ shows the 
effect on the 1-cr confidence level of using different values of 
T. The dashed line represents the 1-er contour computed for the 
correct value of V used for the simulations (T = 0.21); the 
upper and lower contours are computed for V = 0.23 and T = 
0.19, respectively. From this plot, we can see that mistaking T 
by 10% will introduce a systematic of about 16% in Br and 
10% in n . 

5. Conclusion and prospects 

Background quasars and foreground galaxies are cross- 
correlated on angular scales between ~ 1' and ~ 1° through 
the magnification bias due to weak lensing by large-scale struc- 
tures. The effect probes the dark matter distribution in the 
Universe in much the same way as measurements of cosmic 



shear do, but by counting objects instead of measuring shapes, 
and thus without the crucial assumption that the ellipticities of 
background galaxies are uncorrected. In addition, the lensing- 
induced QSO-galaxy correlations probe the relation between 
the distributions of galaxies and dark matter. 

The existence of statistically significant correlations be- 
tween foreground galaxies and background quasars has firmly 
been established. Up to now, however, the quasar-galaxy cor- 
relation function wqg could not be reliably measured be- 
cause homogeneous, sufficiently deep, wide-field galaxy sur- 
veys were missing. With the advent of such surveys, for which 
the Sloan Digital Sky Survey (SDSS) is the prototypical ex- 
ample, a measurement of this weak lensing effect with high 
signal-to-noise ratio is now within reach. 

Using quasar-galaxy cross-correlations and assuming a lin- 
ear and deterministic bias, Bemtez & Sanz (2000) showed 
how to measure ilo/b at a given angular scale from the ratio 
W QG (#) / wgg (0) between the quasar-galaxy cross-correlation 
and the galaxy-galaxy autocorrelation functions. However, 
their method discards the information contained in the ampli- 
tude and angular shape of wqg- In this paper, we investigated 
what kind of cosmological information can be expected from 
an accurate measurement of wqg using the SDSS data. Our 
main assumptions were that 

- galaxies are biased relative to the dark matter, possibly in a 
non-linear or stochastic way. 

- the shape parameter of the power spectrum either satisfies 
T oc f2o h, or can be fixed by the peak location in galaxy 
power spectra; 

- the power spectrum is normalised such that the local abun- 
dance of rich galaxy clusters is reproduced; and 

- the lensing magnification is weak. \fi— 1| <C 1. 
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Under these assumptions, wqg depends in principle on 
the density parameter ilo, the cosmological constant the 
Hubble constant h (or T), the normalisation of the power spec- 
trum <rg, and the bias parameter b. This parameter range is 
too large to allow any accurate cosmological constraints from 
the quasar-galaxy correlations alone. However, the parameter 
range can safely be reduced. We showed that: 

1. the dependence of the quasar-galaxy cross-correlation 
function on the cosmological constant fl\ is entirely negli- 
gible on scales above one arc minute; 

2. assuming the relation r = ilo h, we find an interesting an- 
gular range within which the amplitude of wqg is insensi- 
tive to f2rj> an d where its amplitude is simply proportional 
to 6(40') of; 

3. neglecting any relation between T and the other parameters 
and using the cluster normalisation, the matter density and 
the bias parameter can be constrained with with 17% and 
10% relative accuracy, respectively; 

4. knowing the value of the bias parameter, it is possible to 
accurately measure of Hq and erg, namely with 17% and 
7% relative accuracy, respectively. 

We used a scale-independent bias, but models or measured 
constraints on b(6) can be included into the calculation of the 
cross-correlation function wqg {&)■ 

Our knowledge of the cosmological parameters will greatly 
improve over the next few years, partly through the wide- 
field galaxy surveys themselves, but most dramatically through 
the CMB experiments. Therefore, the quasar-galaxy cross- 
correlations will be less important for determinations of cosmo- 
logical parameters, although they will allow important consis- 
tency checks. However, the most important contribution from 
a measurement of wqg will be a direct determination of the 
galaxy bias factor on a broad range of scales, for different lu- 
minosities and morphological types of galaxies. 

Williams & Irwin (1998) and Norman & Williams (2000) 
cross-correlated APM galaxies with LBQS and 1-Jansky 
quasars and claimed significant galaxy overdensities around 
quasars on angular scales of order one degree. As the ampli- 
tude of their correlations cannot be explained in terms of weak 
lensing by large-scale structures in a CDM universe, it is highly 
important to reexamine them. 

Finally, the shape of wqg depends on the shape of the un- 
derlying dark-matter power spectrum, and a measurement of 
wqg will therefore test whether the CDM power spectrum ade- 
quately describes the weak-lensing effects by large-scale struc- 
ture. This, of course, is also probed by cosmic shear, but only 
under the assumption that intrinsic galaxy ellipticities are un- 
correlated. Measuring wqg thus provides a useful, indepen- 
dent, weak-lensing based method for directly constraining the 
dark-matter power spectrum. 

Acknowledgements. We would like to thank Yannick Mellier, Ludovic 
Van Waerbeke and Simon White for useful discussions. 

References 

Bacon, D., Refregier, A., Ellis, R. S. 2000, MNRAS 318, 625 



Bardeen, J. M., Bond, J. R., Kaiser, N.; Szalay, A. S., 1986, ApJ 304, 
15 

Bartelmann, M., 1995, A&A, 298, 661 

Bartelmann, M. & Schneider, P., 1993, A&A, 271, 421 

Bartelmann, M. & Schneider, P., 1994, A&A, 284, 1 

Bartelmann, M., Schneider, P., & Hasinger, G., 1994, A&A, 290, 399 

Bartelmann, M. & Schneider, P., 2001, Physics Reports 340, 291 

Bemtez, N., Martinez-Gonzalez, E., 1995, ApJ 448, L89 

Bemtez, N., Sanz, 2000, ApJ 525, LI 

Bemtez, N., Sanz, J. L., Martinez-Gonzalez, E., 2001, MNRAS 320, 
241 

Bernardeau, R, Van Waerbeke, L., Mellier, Y, 1997, A&A 322, 1 
Dekel, A. & Lahav, O., 1999, ApJ, 520, 24 
Dolag, K. & Bartelmann, M., 1997, MNRAS 291, 446 
Dodelson, S., Narayanan V. K., Tegmark, M., Scranton, R., Budavari, 
T., Connolly, A., Csabai, I., Eisenstein, D., Frieman, J. A., Gunn, 
J. E., Hui, L., Jain, B., Johnston, D., Kent, S., Loveday, J., Nichol, 
R. C, O'Connell, L., Scoccimarro, R., Sheth, R. K, Stebbins, A., 
Strauss, M. A., Szalay, A. S., Szapudi, I., Vogeley, M. S., Zehavi, 
I., et al., for the SDSS Collaboration, astro/ph 0107421 
Eke, V. R., Cole, S., Frenk, C. S., 1996, MNRAS 282, 263 
Fugmann, W., 1990, A&A, 240, 11 

Guimaraes, A., C, C, van de Bruck, C, Brandenberger, R., H., 2001, 

MNRAS 325, 278 
Freedman, W. L., Madore, B. R, Gibson, B. K, Ferrarese, L., Kelson, 
D. D., Sakai, S., Mould, J. R., Kennicutt, R. C. Jr., Ford, H. C, 
Graham, J. A., Huchra, J. P., Hughes, S. M. G., Illingworth, G. D., 
Macri, L. M., Stetson, P. B„ 2001, ApJ 553, 47 
Heavens, A., Refregier, A., Heymans, C. 2000, MNRAS 319, 649 

Hoekstra, H, Franx, M., Kuijken, K., 2000, ApJ 532, 88 

Hoekstra, H, Yee, H.K.C, Gladders, M., 2001a, |astro-ph/0 107413 
Hoekstra, H, et al., 2001b, in preparation 



Kaiser, N., Wilson, G., Luppino, G. A., 2000, |astro-ph/0003338| 
Kauffmann, G., Colberg, J. M., Diaferio, A., White, S. D. M., 1999, 

MNRAS 307, 529 
Norman, D. J. & Impey, C. D., 1999, AJ 118, 613 
Norman, D. J. & Impey, C. D., 2001, ApJ 552, 473 
Norman, D. J., Williams, L. L. R. 2000, AJ 1 19, 2060 
Peacock, J.A, Dodds, S.J., 1996, MNRAS 280, L19 
Pei, Y. C, 1995, ApJ 438, 623 

Rodrigues-Williams, L. L. & Hogan, C. J., 1994, AJ, 107, 451 
Sanz, J. L., Martinez-Gonzalez, E., Bem'tez, N., 1997, MNRAS 291, 
418 

Schneider, P. 1998, ApJ 498, 43 

Scranton, R., Johnston, D., Dodelson, S., Frieman, J. A., Connolly, A., 
Eisenstein, D. J., Gunn, J. E., Hui, L., Jain, B., Kent, S., Loveday, 
J., Narayanan, V, Nichol, R. C, O'Connell, L., Scoccimarro, 
R., Sheth, R. K., Stebbins, A., Strauss, M. A., Szalay, A. S., 
Szapudi, I., Tegmark M., Vogeley, M Zehavi, I. & the SDSS 
Collaboration, 2001, |astro-ph/0 1074 16 

Sugiyama, N., 1995, ApJS 100, 281 

van Waerbeke, L. 1998, A&A 334, 1 

van Waerbeke, L., Bernardeau, R, Mellier, Y, 1999, A&A 342, 15 
van Waerbeke, L., Mellier, Y, Radovich, M., Bertin, E., Dantel-Fort, 
M., McCracken, H. J., Le Fevre, O., Foucaud, S., Cuillandre, J.- 
C, Erben, T., Jain, B., Schneider, P., Bernardeau, R, Fort, B., 
2001, A&A 374, 757 

Williams, L. L. R., Irwin, M. 1998, MNRAS 29 8, 378 

Wang, X., Tegmark, M., Zaldarriaga, M., 2001, |astro-ph/0 105091 



Wittman, D. M., Tyson, J. A., Kirkman, D., Dell' Antonio, I., Berstein, 

G. 2000, Nature 405, 143 
York, D. G., Adelman, J., Anderson, J. E., Anderson, S. R, et al. 2000, 

AJ, 120, 1607 



